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Preview-Based Stable- Inversion for Output Tracking 


Qingze Zou * Santosh Devasia t 


Abstract 

Stable Inversion techniques can be used to achieve high-accuracy output tracking. 
However, for nonminimum phase systems, the inverse is noncausal - hence the inverse 
has to be precomputed using a pre-specified desired-output trajectory. This requirement 
for pre-specification of the desired output restricts the use of inversion-based approaches 
to trajectory planning problems (for nonminimum phase systems). In the present article, 
it is shown that preview information of the desired output can be used to achieve online 
inversion-based output tracking of linear systems. The amount of preview-time needed is 
quantified in terms of the tracking error and the internal dynamics of the system (zeros 
of the system). The methodology is applied to the online output tracking of a flexible 
structure and experimental results are presented. 


‘Graduate Student, Mechanical Engineering Department, University of Utah, SLC, UT, 84112. 
Email: qingze@eng.utah.edu 

f Assistant Professor, Mechanical Engineering Department, University of Utah, SLC, UT, 84112. 
Email: santosh@eng.utah.edu 



1 Introduction 


Inversion of system dynamics can be used to find inputs which achieve exact output-tracking 
[10, 22]. For systems with nonminimum phase dynamics, inputs found through standard inver- 
sion techniques tend to be unbounded and therefore, cannot to be used for practical output- 
tracking. Recently, stable inversion approaches [5, 11, 1] have been developed for nonminimum 
phase systems, which yield bounded inputs for exact output-tracking. The application of such 
inversion-based output tracking has been studied for several nonminimum systems, like flexible 
manipulators [20, 16], aircraft systems [18, 24, 4], and high precision positioning of piezo-probes 
for nano-technology [2] . A critical difficulty with these inversion-based approaches is that the 
inverse is noncausal (for non- mini mum phase systems) and therefore the desired output tra- 
jectory must be pre-specified. This requirement for pre-specification of the output-trajectory 
can be a substantial limitation on the use of the non-causal inversion-based approach, and it 
limits the inversion-based approach to trajectory planning applications. The main contribution 
of this paper is to show that, for linear systems, the noncausal inverse can be computed using 
a preview-based approach, which allows the inversion process to be applied online. This article 
also quantifies the amount of preview-time needed in terms of the desired tracking accuracy 
and the location of the zeros of the system. Implementation issues are discussed and the ap- 
proach is experimentally verified by applying it to the output tracking of a flexible structure 
with nonminimum phase dynamics. 

A major result in output tracking is the solution of the output regulation problem for 
linear systems by Francis [9]. These results were generalized for the nonlinear case by Isidori 
and Byrnes [13]. The desired outputs are assumed to be generated by an exosystem and the 
linear regulator is easily designed by solving a manageable set of linear equations. A problem, 
however, with the regulator approach is that the exosystem states are often switched to describe 
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the desired output; this leads to transient tracking errors after the switching instants. Such 
switching- caused transient errors can be avoided by using inversion-based approaches to output 
tracking [5, 7]. Thus, it is advantageous to use inversion-based output tracking when exact- 
tracking of a particular output trajectory is required. In the inversion based approach, the 
system dynamics [10, 22] is inverted to find the input that exactly tracks a single specified output 
trajectory (rather than track a class of outputs as in the case of the output-regulator). Inversion 
for nonminimum phase systems is challenging since the standard approaches lead to unbounded 
inputs [10]. Stable inversion techniques resolve this problem of unbounded inverse-inputs by 
finding bounded (but possibly noncausal) exact-tracking input-state trajectories [5, 11], The 
noncausality of the inverse implies that the entire output trajectory needs to be known ahead 
of time which restricts the use of inversion-based approaches (for nonminimum phase systems). 
This motivates the current work, which shows that the noncausal inverse can be found by 
using preview information and thereby enables the online implementation of the inversion- 
based output tracking technique for nonminimum phase systems. 

Other approaches have also used preview information of the desired output trajectory for 
output tracking, for example, to alleviate the problems due to nonminimum phase dynamics 
by exploiting actuator redundancy [25]. Another use of previewed information of the output 
is in linear quadratic-based (lq-based) optimal output tracking controllers (see [17], Chapter 
4). Tomizuka and coworkers (e.g., [23]) have shown that the performance of finite-time-preview 
controllers approach the performance of the infinite-preview controller as the amount of preview 
time increases. In these works, the goal is to trade-off the tracking requirement to reduce the 
magnitude of the inputs. In contrast, inversion based approaches aim to achieve high accuracy 
control of the desired output trajectory. Trading off the accuracy output-tracking requirement 
to achieve other goals like vibration minimization is also possible within the inversion-based 
framework (see [8]). The resulting inverse controller is noncausal, which can also be imple- 



mented using the preview-based controller discussed in the current article. 

The paper is organized in the following format. The inversion-based output tracking scheme 
and its dependence on the solution of the system’s internal dynainics is presented in Section 
2. The inversion problem is then solved using preview information of the output, and imple- 
mentation issues are studied in Section 3. In section 4, the preview-based inversion approach 
is applied to the output tracking control of a flexible structure and experimental results are 
presented. Discussions are in section 5 and our conclusions are in Section 6. 


2 Stable Inversion for Nonminimum Phase Systems 

In this Section, the inversion-based output tracking approach is presented. It is shown that 
solving the inversion problem is equivalent to finding bounded solutions to the systems ’s internal 
dynamics. 


2.1 Output tracking using inversion of system dynamics 

Consider a linear system described by 


x(t) = A x(t) + B u(t) 
y(t) = C x(t) 


which has the same number of inputs as outputs, u(t),y(t ) E 3£ p , and x(t) E 3? n . We assume 
that the system is stabilizable. Let y d ( ) be the desired output trajectory to be tracked. Then 
in the inversion-based approach we, first, find a nominal input-state trajectory, [u ff (-),x ref (-)] 
that satisfies the system equations (1) and yields the desired output exactly, i.e., 


•fVe/(0 
Vd(t ) 


■> 

A x Tef (t) + B u ff (t) 


V£ E ( — oo, oo) 


( 2 ) 
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and, second, we stabilize the exact-output yielding state trajectory, x re f (• ) , by using state feed- 
back (see Figure 1). Thus x{t ) — > x re j{t ) and y(t) — * y d {t) as t — * oo and output tracking is 
achieved. While stabilization of the reference state trajectory can be easily achieved through 
standard techniques [14] like state feedback of the form K[x(t) — x re f(t)], the main challenge is 
to find the inverse input-state trajectory [w//( - ), x re f(-)] — especially for systems with nonmini- 
mum phase dynamics. This paper addresses the on-line computation of the inverse input-state 
trajectory using preview information of the desired trajectory, y d . 


Feedforward 



Figure 1. The output tracking control scheme. 

2.2 Stable inversion scheme 

In this subsection, it is shown that finding the inverse input-state trajectory is equivalent to 
finding bounded solutions to the system’s internal dynamics. Let the linear system (1) have 
a well defined vector relative degree, r := [jq, r 2 , . . . , r p ]. Then the output’s derivatives are 

given as: 

d rk vt 

-Jt 7T = C ^ Arkx + C k A r *~ x Bu (3) 
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where C k is the k th row of C, and 1 < k < p. In vector notation let equation (3) be rewritten 
as 


y (r) {t) = A x x(t) + B y u(t) 


( 4 ) 


where 


y 


(r) 


d ri Vi d r2 y 2 
dt r ' dt r 2 


d rp y P 

dt T p 


1 T 



CiA l ri 


CxA^-'B 

4t:= 

C 2 A l" 2 

; By 

C 2 A r ^B 


C p A r v 


C p A r p~ l B 


and By is invertible because of the well-defined relative degree assumption. Equation (4) mo- 
tivates the choice of the control law of the form 


u ff (t) = By 1 [yd (r) (0 _ A x x(t)] 


( 5 ) 


for all t € (— oo, oo). Substituting this control law in equation (4), it is seen that exact tracking 
is maintained, i.e., 

y (T> (t) = vP (()• 

To study the effect of this control law consider a change of coordinates T such that 


m 

v{t) 


= T x(t) 


(6) 


where £(£) consists of the output and its time-derivatives 

dtr 2 -i 


r . d ri ~ 1 y l . d r *~ l y 2 

at) ■= [yu 2/1, ■ • • , -j—r, 2/2, z/2, ■ • - , 


dr>- 


<r*- l vp v 

’ dt T p~ l J 
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The system equation (1) can then be re-written in the new-coordinates as 


€(t) — + A 2 T) + B\U 

fj{t) = A 3 ^ + A4T] + B 2 U 


(7) 

(8) 


where 


1 

Ax 

A 2 

; and B := 

Bx 


A 3 

1 

TT 


. & 


= TB 


In the new coordinates, the control law for maintaining exact tracking (Equation (5)) can be 
written as 


u ff (t) = B y 1 [yA r] (t) - A ( £ d (t) - Ar,r)(t)] 


(9) 


where 


* 

A„ 


:= A x T~\ 


s 

was 


Note that the desired £(•) is known when the desired output trajectory y d {-) and the output 
time derivatives are specified. This desired £(■) is defined as &(■). Since the control law 
chosen such that exact tracking is maintained, y {r) {t) = y d \t) we also have £(t) = £ d (t), and 
Equations (7) and (8) become 


«<) = 

r](t) = A 3 £ d (t) + A 4 7)(t) 4- B 2 B~ l [y/ r \t) — A^ d (t) — A v T}(t )] 
: = Afj rj(t) 4 - B n Y d (t) 


(10) 


( 11 ) 


where 


K : = ^4 - fyB^Ar, ; B n := 

B 2 By X A 3 - B 2 B~ 1 aA 

; and Y d := 

2/5(0 


L J 


1 

1 


( 12 ) 
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This is the inverse system, and in particular, Equation (11) represents the internal dynamics. 
If a bounded solution, %(•), to the internal dynamics (11) can be found then the exact-tracking 
feedforward input can be found through equation (9) as 


«//(*) = B; 1 [y d (r) (0 - A&(t) - A v rj d (t )] 


(13) 


and the reference state trajectory can be found as 


x ref {t) = T 1 


Ut) 

Vd(t ) 


(14) 


Thus a bounded solution to the internal dynamics (11) is required to find the inverse for applying 
the output tracking scheme shown in Figure 1 . 


2.3 Bounded solutions to the internal dynamics 

We restrict the following discussion to systems with hyperbolic internal dynamics, i.e., none of 
the zeros of the system (Equation 1) lies on the imaginary axis of the complex plane (a technique 
to address the nonhyperbolic case can be found in [4]). If the internal dynamics is hyperbolic, 
there exists a transformation U such that the internal dynamics (11) can be decoupled into a 
stable subsystem (< 7 S ) and an unstable subsystem (< 7 U ): 


d a {t) = A a a s {t ) + B s Y d (t ) (15) 

0u(O = A u a u (t) + B u Y d {t) (16) 


where 


a{t) 


° s {t) 

Gu{t) 


= u v {t) 


(17) 


Bounded solution to the internal dynamics in the transformed coordinates can then be found 
as 
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<r a (t) = f e A ‘^B s Y d {r)di 

J -oo 


( 18 ) 


/ oo 

e- A ^B u Y d (r)d T (19) 

In the new coordinates, the feedforward control law in Equation (13) can be written as: 

«//(*) = B - 1 [y r d (t) - AsUt) ~ A v U s a s (t ) - A n U u a u {t)\ (20) 

where U~ 1 := \(j s U u is partitioned according to the partition of a in Equation (17). 

This completes the inversion technique. To summerize: the bounded solutions found through 
Equations (18) and (19) are used to find a bounded solution to the internal dynamics, r) d , by 
using Equation (17). The inversion is then completed by finding the reference state and input 
trajectories by using Equations (14) and (20), which are then used in the control scheme shown 
in Fig. 1 to obtain output tracking. Note that only the past information is needed to compute 
the solution to the stable subsystem of the internal dynamics by using Equation (18). However, 
to find a bounded solution to the unstable subsystem ( cr u ) by using Equation (19), the desired 
output must be completely specified (including future information). This is the main problem, 
which restricts the use of inversion to trajectory planning (where it may be acceptable to solve 
Equation (19) off-line). 


3 Preview Based Inversion 

In this section, we discuss the online computation and implementation of the inverse using 
preview-information of the desired output trajectory, which enables the tracking of output 
trajectories that are specified on-line. We begin by quantifying the relationship between the 
preview-time and tracking error. 
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Let the desired output, y d (and its time-derivatives) be given for a preview-time of T p 
seconds, i.e., at time t, Y d (r) (defined in Equation (12)) is known for all t < r < t + T p . This 
preview information is used to approximate the solution to the internal dynamics (in particular, 
to approximate the bounded solution to the unstable subsystem, cr u , given by Equation 19). 
The approximated solution, a u is found as 

rt+Tp 

**(t) = - j e~ A ^B u Y d (r)dr ( 21 ) 

Let the error between the exact-solution to the unstable subsystem found through Equation 
(19) and the approximated- solution found through Equation (21) be defined as: 

- d u (t) (22) 

We show in the following Lemma that e Cu (t) can be made arbitrarily small by having a large 
enough preview time. Furthermore, we show that the error in computation of the inverse input 
by using the finite preview also converges to zero as the preview time, T p , increases. The error 
in computing the inverse-input, e u (t), is defined as 

e u (t) = u ff (t) - u ff (t) (23) 

where U//(t) denotes the input obtained with infinite preview using Equation (20) , and Uff{t ) 
denotes the preview input obtained with finite preview when a u is replaced by the approximate 
solution < 7 U in Equation (20) 

u/f(t) := B~ l |yj(£) - A^ d (t) - A v U s a s (t) - A v U u d u (t)j . (24) 

This implies that, as the preview time increases, the approximated input trajectory approaches 
(arbitrarily closely) the exact output- tracking inverse trajectory. 

Lemma X Let the desired, trajectory and its derivatives be bounded , i.e. there exists a positive 
scalar Me R, such that ||Y d (t)|| 2 < M for all time t. Then for any scalar e > 0, there 
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exists a finite time T* such that the error ( in computing the inverse input by using the preview 
information of the output) can be made smaller than c if the preview time is larger than T ) , i.e., 
||cu(^)ll 2 ^ £• Itn the above notation, \\z \\2 is defined as the standard Euclidean norm for any 
vector z. 

Proof: Since the subsystem (16) is unstable, — A u is Hurwitz, and therefore positive scalars a 
and k can be chosen such that [3] 

\\e~ Au % < Ke ~ at V {>0 (25) 

where, given a matrix F € R n * n , ||.F|| 2 denotes the induced matrix norm, defined as 

||F|| 2 := sup ||Fx|| 2 . From the definition of e Cu in Equation (22) 

M2=l 


ll^Wlla = 


< 

< 

< 


<Tu(t) - d u (t)|| 2 


roc 

/ e _ ^“ (T "‘ ) B u Y d (r) dr|| 2 (using Equations (19) and (21)) 

« f+Tp 
oo 




2 ||£ u Y d (r)|| 2 dr 


roc 

J t+T p 

roc 

Ki / ||e^ (T "°|| 2 dr (with K x := M\\B U \\ 2 ) 

Jt+T p 

~ e~ aTp (using Equation (25)) 


(26) 


From the definition of the computational-error in finding the inverse-input (Equation (23)) and 
from Equations (20), (22) and (24), we have: 


sup ||e u (f)|| 2 


sup II B~% \\A n \\ 2 \\U u h ll^ u (t)|| 2 


:= supJf,|K(«)lla (*a := ||-B“ x || a N^llall^lls) 

1C 


< — e 

a 


OcTrt 


(using inequality (26) and fC := KiK 2 k) 


(27) 



For any given e > 0, choosing 


t;> 


In rr 


a 


(28) 


and substituting any T p > T* in Equation (27) completes the proof. □ 

It is noted that for proving this lemma, it is sufficient for ||.B u Y (f (-)|| 2 to be bounded, which can 
be less restrictive than requiring ||Y d (-)|| 2 to be bounded. Using similar arguments, it is also 
possible to show that the error in computing the inverse reference-state-trajectory, x re f, can be 
made arbitrarily small by having sufficiently large preview time. 

The next lemma shows that the error in output-tracking, due to errors in computing the 
inverse, can also be made small by choosing a sufficiently large preview time (for a similar 
argument, see [6]) 


Lemma 2 Let the original system (1) be stable (or stabilized with feedback before the inversion 
is applied). Then the tracking error, on applying the finite-preview- based input (u ff ) as a 
feedforward input, can be made arbitrarily small by choosing a large enough preview time. 

Proof: From Eq. (1), the dynamics of the state-error e x (t) := x(t) - x ref (t), where x(t) is the 
system-state when the finite preview input u ff {t) is applied to system (1), can be described by 

&x{t) — A e x (t ) + B e u (t) (29) 

The error in output-tracking e y (t) := y(t) - y d (t) = C e x (t) can then be bounded as 

IMOIIa = IIC^WII* < ||<?||a IM*)lla 

^ ll^lb J ||e i4(t_r) || 2 ||B|| 2 ||e u (T)|| 2 dr (using solution of Equation (29) ) 

J— OC 

< l|C|| 2 ||B|| 2 ( S up||e„(T)|| 2 ) /' ||e' 1( ‘- r) || 2 dr (30) 

If the original system(Ecjiiation(l) is stable (or if the inversion is carried out after stabilizing the 
system), then there exit real positive numbers M\ and /? such that ||e y4 ( t_T ^|| 2 < Mi e -/3 ( £ ~ T ). 
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Substituting this in Equation (30) yields 








IMOIk < ||C|| 2 ||5|| 2 (sup ||e u (r)|| 2 ) f M x e ^ T) dr 

T J — oo 

Mi 

= M 2 (— ) (sup ||e u (r)|| 2 ) ( with M 2 := ||C|| 2 ||B|| 2 ) 

MiM 2 A: _ c 


C — ‘--^- e -aTp 


Q/? 


(using inequality (27)) 


( 31 ) 


The tracking error, ||e I/ (f)|| 2) can be kept smaller than e max for all time, t, by choosing the 
preview time, T p , as 

In m,m 2 k 

T p > T; > - Q ^ maT (32) 

which completes the proof □ 


Remark: The premew time T p needed to achieve a desired accuracy in output tracking depends 
on the distance of the right half plane zeros of the system (Equation (1)) from the imaginary 
axis of the complex plane. As this distance, increases, a in Equation (32) can be chosen larger 
(defined in Equation (25)) and therefore, a smaller preview-time can be chosen to achieve the 
desired accuracy in output-tracking. 


4 Example: Flexible Structure Control 

v 

The experimental flexible structure considered here consists of two discs which are connected 
by a thin freely rotating shaft as shown in Figure (2). The system input, U(t), is the voltage 
(Volts) applied to a DC motor, and the output, 0 2 , is the angular rotation (in degrees) of the 

V ' 

disc which is farther from the motor. The system equations (obtained experimentally with a 
HP3562A Dynamic Signal Analyzer) can be written in the following state-space form where the 
state vector, x is chosen as 

X [#1 9y 02 02^ 
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0 

1 

0 

0 


0 


-3.656 

-0.436 

3.573 

-0.091 


21.9027 

X = 

0 

0 

0 

1 

X + 

0 


3.245 

-0.126 

-3.259 

-0.076 


3.588 


u 


y := 


Ax + Bu 
62 = [0 0 1 0] x 



(33) 


(34) 


Figure 2. The experimental flexible structure. 

4.1 System inverse 

The relative degree of the above system is two and hence the output has to be differentiated twice 
to relate the input and the output (as described in Section (2.2) with £ := [y, y] r = [0 2 , fW 7 ')' 
Note that £ is known when preview information of the desired output and its derivatives are 
defined. The internal dynamics of the system are described by 77 := [ 0 j, 6 1 ] T . The inverse input 
Uf f can be written as (using Equation (13) and the last row of Equation (33)) 


1 


«// ( 0 = 3^8 {&(<) - ; [3.245 - 0.126] rj(t ) - [-3.259 -0.076]&(t)} 


(35) 


and the internal dynamics is given by substituting this control law into the first two rows of 
Equation (33) 
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0 

1 


0 

0 

0 

V = 

-23.4658 

0.3306 

V + 

23.4658 

0.3753 

6.1041 


Y<t 

0.3306 J |_ 234658 0.3753 6.1041 J 

:= A^t) + BrjYd 

where, Y d = [y d y d y d ] T and the entire internal dynamics is unstable, i.e., a u = rj 


(36) 

(37) 


4.2 Y d generation 

There are several methods available to generate the desired trajectory y d and its time derivatives. 
One method is to predict the future y d by using polynomial extrapolations of the past desired 
trajectory signals [19] and then differentiating these polynomials. Another approach is to use a 
command generator [23] and then switching its states to generate online changes of the desired 
trajectory. For our experimental system, we assumed that the future desired trajectory y d 
could be obtained for the preview time T p (a potentiometer was adjusted to define the desired 
output). Note that to compute the inversion-based control, the output’s time derivatives must 
also be specified. Although these derivatives, y d , y d , could be found by direct differentiation, we 
avoided noise problems associated with direct differentiation by using a biquad filter [12] shown 
in Figure 3. This filter was chosen with a cut-off frequency of 1 Hz to match the modelled 
bandwidth of the flexible structure - only the first two modes of the flexible structure were 
modeled and the higher-frequency dynamics are not represented by Equation (33) - see [8] for 
details of the flexible structure modeling. 


14 




Figure 3. On-line generation of the desired trajectory and its time derivatives [12]. 

4.3 Online implementation of preview-Based inversion 

The preview-based solution to the unstable subsystem of the internal dynamics, i.e., the inte- 
gral Equation (21), was computed online by discretization by using the fourth-order Simpson 
formula [15] 

&u(t) = - J e- A ^ T ^B u Y d (r)dr 

h N 

« 3 i) (38) 

° i= 0 

where the sampling time is T a such that N := T p /T s is an integer, Y d (t, i) = Y d (t + i * T s ) 
and: 

1 -e~ Au B u if i = 0 

—2 * B u if i > 0 and i is even 

-4 * e~ Au{ - i * T ‘' ) B u if i is odd 

( 39 ) 
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Note that the matrices M{i) can be precomputed and stored to reduce the online computa- 
tions. This computation-scheme is shown in Figure 4. The computation-error, due to time- 
discretization , can be reduced by choosing a small enough sampling time(T s ). 



Figure 4. The schematic of on-line calculation of the unstable part of internal dynamics. 

4.4 Experimental Results 

Experimental results for output-tracking with two-different preview- time, T p = 20 seconds and 
Tp = 50 seconds, are shown in Figure 5, which illustrates the improvement in output-tracking 
as preview-time increases. The experiments include feedback-stabilization, which was added 
to account for unmodeled dynamics in the system (like friction in the experimental system 
and static imbalance of the discs, which created a tendency in the discs to settle in a specific 
orientation). It is noted that trajectories for the two cases are different in Figure 5 because 
these desired output trajectories were generated on-line for the two different preview-time cases. 
As shown in Figure 5, preview-based inversion improves tracking performance with increasing 
preview time, and online specification of the output trajectory is possible. 










Figure 5. Experimental results. The solid line is the desired trajectory, the dashed line 
represents the output trajectory. Plot (a) is for 20 seconds preview time, and plot (b) is for 50 

seconds preview time. 


5 Discussion 

For nonminimum phase systems, recent stable inversion-based approaches can be used to achieve 
high-accuracy output tracking. However, the entire desired output-trajectory has to be pre- 
scribed for the off-line computation of the inverse — this is a significant limitation since the 
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desired trajectory cannot be changed online. This limitation has been alleviated by the current 
preview-based approach which allows the online implementation of the inversion-based output- 
tracking controllers. The methodology allows the application of the inversion-based approach 
to systems like flexible manipulators and servo-positioning systems where it is necessary to 
change the desired trajectory. Such online-changes are allowed by the current technique if a 
preview of the desired trajectory is possible. This requirement for preview information is not a 
drawback of the current controller. Rather, this is necessary for high-accuracy output tracking 
because there are performance limitations on output- tracking for nonminimum phase systems 
if output-preview is not available [21]. An important result of the preview-based inversion ap- 
proach is that the output-tracking error can be made arbitrarily small with sufficiently large 
preview time. The quantification of the amount of preview time in terms of the location of 
system-zeros in the imaginary plane(see section 3) can aid in the initial design of a system. In 
particular, the preview time needed can be made smaller by choosing design-parameters such 
that right-half plane zeros (if any) are far away from the imaginary axis of the complex plane. 
In summary, the preview-based approach will help in the implementation of inversion-based 
control laws for high-accuracy, on-line output-tracking of nonminimum phase systems. 

6 Conclusion 

We have developed and implemented preview-based output tracking using on-line inversion 
for non-minimum phase systems. The preview time needed was quantified in terms of the 
required accuracy in output tracking, and related to the system-zeros. Implementation issues 
were discussed and the technique was verified by applying it to the output tracking of an 
experimental flexible structure. 
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